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)Cl ' Abstract 

. The Davey-Stewartson I equation is a typical integrable equation in 2+1 dimensions. 

I Its Lax system being essentially in 1+1 dimensional form has been found through nonlin- 

earization from 2+1 dimensions to 1+1 dimensions. In the present paper, this essentially 
1+1 dimensional Lax system is further nonlinearized into 1+0 dimensional Hamiltonian 
systems by taking the Bargmann constraints. It is shown that the resulting 1+0 dimen- 
sional Hamiltonian systems are completely integrable in Liouville sense by finding a full 
set of integrals of motion and proving their functional independence. 
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^ . 1 Introduction 

The Davey-Stewartson I (DSI) equation is a famous 2+1 dimensional integrable equation which 
describes the motion of water wave This equation has localized soliton solutions and has 
been studied in various ways, such as inverse scattering 0, |^ , binary Darboux transformation 
0], nonlinearization to 1+1 dimensional problems |^ etc. 

For 1+1 dimensional integrable systems, the nonlinearization procedure, both mono-non- 
linearization and binary nonlinearization f^, reduces them to finite dimensional (1+0 di- 



mensional) integrable Hamiltonian systems H, |T0|, ITT], O]. Therefore, it transforms a partial 



differential equation to a systems of ordinary differential equations. This greatly simplifies the 
procedure of getting solutions, at least numerical solutions. Some important explicit solutions, 
especially periodic or quasi-periodic solutions have been obtained in this way. 

The nonlinear constraint method has also been applied to some 2+1 dimensional equations 



like the KP, MKP, N-wave equations etc. [13, O, n5|. For the DSI equation, we have already 



found its new Lax system ([l|) by nonlinearization in which all the derivatives are separated. 
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Each pair of equations in this system is 1+1 dimensional. Hence the derived Lax system is 
looked as essentially 1+1 dimensional because we can use 1+1 dimensional method to solve 
it. It is possible to nonlinearize this essentially 1+1 dimensional system again to get finite 
dimensional Hamiltonian systems. 

In the present paper, we show that there are Bargmann constraints which reduce the DSI 
equation to finite dimensional Hamiltonian systems. These Hamiltonian systems have a full set 
of integrals of motions and these integrals of motion are functionally independent in a dense 
open subset of the phase space. Therefore, these Hamiltonian systems are completely integrable 
in Liouville sense. 

2 Nonlinearization 



We consider the following Lax system 



^ iX i/ ^ 
iA ig 

^ iA u if 

—u — iA —ig 







(1) 



-2iX'+i\u\''+ivi -2u\+iuy -2i/A-2/; 
2'uA+i-Uj^ 2iA^ — i|Mp — if2 2ig\ — 2gy 
-2if\ + 2fy 2igX + 2gy -2i{\f\^-\g\^)) 
Here u, f and g are complex functions, Vi and V2 are real functions. 



^yx, ^xt 



Its integrability conditions 
three parts. 

(1) DSI equation 

—iut = Uxx + Uyy + 2|-up'u + 2(t>i + V2)u 

Vl,x - = V2,x + V2,y = -{\u\'^)x- 

(2) Standard Lax pair of the DSI equation 



$te and ^yt = ^ty consist of the following 



(2) 



Fx + 




u 



u 



+ 2i 



V 



u 
-u 



(3) 



Ux + U. 



—Ux + Uy —\u\—2v2 
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where F = (/, g)'^. 

(3) Nonlinear constraint 



FF* = - 
2 



(4) 



Hence any solution of (^-(^) gives a solution of the DSI equation. 

Notice that if $ is a vector solution of (0) for real A, then = i$ is a solution of the adjoint 
equations 

where each entry of $ is the complex conjugation of the corresponding entry of $. 

In order to obtain the finite dimensional Hamiltonian systems, we first nonlinearize the 
y-equation of (|I]) in the following way. 

Consider the pair = W^^^ and $t = W^^. Let w = {u, —u, if, if, —ig, —ig) containing all 
the variables in |[10|. Then the recursion relations of this AKNS system can be expressed 
in Lenard form 

JGi = KGi.i (/ = 1,2,---) (6) 

where (J, K) is the Lenard pair (J is a non-degenerate constant matrix) which was given by 
PI (here ai, a2, a^, pi, p2, in are ai = i, a2 = -i, as = 0, pi = -2i, P2 = 2i, p^ = 0) 
and {Go, Gi, G2, ■ ■ ■} is the Lenard sequence. The first element of this Lenard sequence is given 
by as 

Go = -2{-u,u,if,if,-ig,-ig) (7) 

which is a kernel of K. 

On the other hand, the variation of spectral parameter can be computed by the general 



formula ITT 



— =Cotr U^^^— 
ow \ ow J 

= Go{i4>l4>2, 10201, 10103, 10301, 10203, 10302 



(8) 



where Co is a constant. 

Now let Ai, ■ ■ ■ , Aat be distinct non-zero real numbers, (0ia, 02a, 03«)'^ be the correspond- 
ing solution of the Lax system (P for A = A^- Let A = diag(Ai, ■ ■ ■ , Xn), = (0ji, ■ ■ ■ , 4>jN 

By the general formulation of nonlinearization, we impose the nonlinear constraint 



T 



which gives the relations 

(<l>2,$i)=-in ($3,$i) = / ($3,$2) = -^? (10) 

where (Vi,V2) = V{V2 for any two vectors Vi and V2. These are the Bargmann constraints 
between {u,f,g) and ($1, $2, 'J's)- 



3 



Remark 1 If we consider the pair = W^^ and = Gq is different. In that case, we 

can not obtain Bargmann constraints, but only Neumann constraints j^^. 



Let 



L(A) 



N 



+ E 



a=l ^ 







V 



'laVla <i>2a<i>la <i>'ia<i>la 
hQ:02a 02Q02a 03q:02q 
Ha^Sa 4>2a4''ia 4'3a4'3a 



Lemma 1 The Lax equations 

= [W^, L] Ly = [Wy, L] U = [W\ L] 



hold if and only if ( tl^ j holds. 

Proof. Let = {4>ia, 4>2a, (p3aV, then 



N 



Since Fa,y = Wy{K)F^ and (Vr^(A^))* = -Vr^(A«), we have 



A- A, 



-[wy{x^),F^F: 



[wyix),Li\)-c] 



-E 



A — Aq, 

[wy{x),L{x)-c]-i 



[wy{x)-wyiXo.),F^F^ 



Hence Ly{X) = [W^^(A), L(A)] if and only if 



[c,wyix)]. 



;i2) 



(13) 



(14) 



(15) 



Written in the components, it becomes (10). If ( p!OD holds, the other two equations in ([T2|) are 
obtained similarly. The lemma is proved. 
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From the above constraints, we have 

= i($3, A$i) + i/(($3, $3) - (<^'i, <^'i)) - ug 
= -i($3, A<l>2) - i^(($3, '^3) - {'^2, $2)) - uf 
fy = i(<l>3, A$i) + i/((<f 3, $3) - ($1, $1)) 

gy = i($3, A<l>2) + i^(($3, $3) - ($2, $2)) (16) 
Ux = '2fg 

Uy = -2($2, A<l>i) + ($2, $ i)((<f 1, $1) - ($2, $2)) 

Vi = 2($i,<l>3)(<l>3,$i) V2 = 2(<l>2,$3)($3,'^'2) 

With the constraints (10), the system (|l]) is changed to a system of ordinary differential 
equations 

$1,^ = iA$i + i/$3 <^2,x = iA$2 + i^$3 

^l^y = iA^i +M<I>2 + i/$3 

^2,y = -u^i - iA<l>2 - ig^3 
$3,^ = i/$i - i^$2 

$1,^ = (-2iA2 + i|M|2 + iv,)^, + (-2mA + my)^2 (17) 
+(-2i/A - 2/,)<|.3 

$2,4 = (2mA + iUy)^! + (2iA2 - i|M|2 - W2)^2 

+ {2igA-2gy)<!>s 

$3,t = (-2i/A + 2/,)<l>i + (2i^A + 2^,)<1>2 

-2i(|/p-|^P)$3 

where u, Uy, vi, V2, f, g, fy, gy are given by (p!oD and (|l^). 

Re0ja and Im^ja (j = 1,2,3; o; = 1,---,A^) form a system of coordinates of the phase 
space R^^. For simphcity, we use the complex coordinates (pja and (j)ja instead of Re(f)ja and 
lm(f)ja- R^^ has the standard symplectic form 

3 N 

^ = 2 I] 51 d lm(0ja) A d Re{(j)ja) 

j=l a=l 

3 N _ 

= i^^I^d^ja Ad^ja. (18) 
The corresponding Poisson bracket of two functions C, and rj is 

By direct computation, we have 
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Lemma 2 ( |7^ is equivalent to three Hamiltonian equations 



50ia <90ja 



where 



H-= -($i,A<l>i)-(<l>2,A<f2)-K$l,$3) 



|2 



+ |($2,$3)P 
i/^ = -($1, A<l>i) + ($2, A$2) - |($1,$3)P 

-|($2,$3)P-K$1,$2)P 
if*= 2($i,A2$i) -2(<I>2,A2$2) 

+4 Re(($i,$3)($3,A$i)) 

+4Re(($2,A)<l>3)(<f3,A$2)) 

+4Re((<l>i,$2)($2,A$i)) 

+2((<l>3,<f3)-(<fl,$l))|($l,$3)P 
-2(($3,<f3)-($2,<f2))|($2,$3)|' 
-((<l>l,$l)-($2,$2))|(<fl,$2)p. 

3 Integrability 

Lemma 3 For any two complex numbers \, fi and any positive integers k, I, 

{trL^(A), trL'(/i)} = 0. 

Proof. Denote e,- = (0, ■ ■ ■ , 0, 1, 0, ■ ■ ■ , 0)^, then 

j 

1 dL{X)_ 1 

A — A„ ^ d(j)ja A — Aq " 



Since 



tr 



A; tr 



V 0(Pja ) \ <JHJja 
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we have 

-^{trL^(A),trL'(/.)} 



kl 



a,b,j,a ^c, - Aa 

- E ^^(^'-'(/^))a.^aayV(^'''(A)W.a (25) 
a,6j,a - Aa A - 

a,b,a A - A« /X - A« 

= E 7— t(t— IT - -7— irl^-'/'&a- 



/i — A V A — fi — Xq^ / 
^J^^^ (^^^^^ " L(/.))[L^-^(A),L^-i(/.)]) = 0. 

The lemma is proved. 

Using this lemma, we can construct involutive integrals of motion from trL'^ {k — 1, 2, • • •). 
For any complex number ^, let 

det(e - L{\)) = e - Piwe + - (26) 

then Pfc(A) is the sum of all the determinants of the principal submatrices of L[X) of order k. 
Suppose the eigenvalues of L{X) are i'i{X) , i'2{X) , ^^{X) , then 

trL^(A) = E^^(A) 

^■=1 (27) 

Pk{x)^ E ^hW---^3kW- 

l<ji<--<jfc<3 

Since Pk{X) can be expressed as a polynomial of tr L'(A) (Z = 1, 2, • • •), {pj(A),pfc(/i)} = for 
any two complex numbers A and /i and two integers j, k > 0. 
Expand Pk{X) as a Laurant series of A 

CXD 

Pk{X) = E E^Jl^X—' (28) 

m=— 1 
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which is convergent when |A| > max |Aq|, then 

l<a<Af 



E'^) = -($i,A'"$i) + ($2,A'"$2) 



+ E E 

i<«<i<3 1=1 



($„A'-i$,) ($,,A™-'%) 



£;(,3) = -($3,a™$3) 

($i,A^-i$i) ($3,A'"-'$i) 
($i,A'-i$3) ($3, A— '$3) 

($2,A'-'$2) ($3, A— '$2) 

(<I>2,A'-1$3) ($3, A— '$3) 

($1, A*$i)($2, A-'^iX^g, A'^$i) 
($1, A''$2)($2, A-^$2)($3, A'^$2) 
($1, A^$3X*2, A^^3X*3, A'=*3) 



-E 

1=1 

m 

+E 

1=1 



+ E 

i+j+fe=m— 2 
ij,fe>0 



(29) 



The sums are zero if the lower bound is greater than the upper bound. These E^)'s are m 
involution. 



Let 



Q2 

^^3 



then 



($i,$i) 
($2,^2) 
(^3, ^3) 

.(1) 



(4^^ - e'^' + 4^V2 

(^?) + + £;?V2 

,(3 




(30) 



-E['> - e[^^ - (Qi - Q2)^^3 



+(01 + 02 - 2^3)//^ + (Ol - 02)^^ 

which have all been expressed as polynomials of -E^^'s. 
Hence the following lemma holds. 



(1) 



(31) 



Lemma 4 



and 



(32) 
(33) 

(34) 



for all k,l — 1, 2, 3; m,p — 0, 1, 2, 
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Now we consider the independence of E^^^s. 
Lemma 5 E^f^ {1 < k < 3; < m < N — 1) are functionally independent in a dense open 



subset ofR^^. 



Proof. Let Pq in R be given by (j)ja = e (j = 1, 2, 3; a = 1, 2, ■ ■ ■ , iV) where e is a small real 
constant. Then, at Pq, 

BES'^'i 

-Ji^ = eX^bk, + Oie') (35) 



where 

/ 1 1 1 \ 

-1 1 

V -ly 



The Jacobian determinant 
J 



PKwi^) p^(l) p(2) p(2) p(3) p{3) ^ 



^l<Q</3<Ar 

J 7^ at Pq if £ 7^ is small enough. Since J is a real analytic function of 



,6N 



J 7^ in a dense open subset of R . Therefore, the Jacobian determinant of 

r(3) p(3) N 



(J^W p(2) p(2) p(3) p(3) 



to 6A^ real coordinates 

( Re(0ii), ■ ■ ■ , Re(03iv), lm(0ii), ■ ■ ■ , Im(03Ar)) 

is of full rank 3A^. The lemma is proved. 
In summary, from Lemma 2-5, we have 



(36) 



9(011, ■ ■ ■ , 4>1N, 021, ■ ■ ■ , 027V, 031, ■ ■ ■ , 037V) 

/ \ 3 (37) 

2^''' n (V-A„) +0(e=^^+2). 



Theorem 1 The Hamiltonian systems given by the Hamiltonians are involutive and com- 
pletely integrable in Liouville sense. The integrals of motion are E^"* [k = 1,2,3; m = 
0, ■ ■ ■ , — 1) given by ^^), which are in involution and functionally independent in a dense 
open subset of the phase space R^^. Moreover, each solution of these Hamiltonian systems 
gives a solution of the DSI equation. 
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